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ABSTRACT 

An approaCli to orbital aechanics, wnicn i5 accessible 
t:: beginnina physics students and presupposes no jcnowieuje of 
calcalus, is presented* A theory of orbits is deve^opea lor tne 
i^verse-sguare central force law which dif f eis* consiueraoiy froii the 
usual deductive approach. This dociiiD^nt begins with gualxtative 
aspects Qt solutions, and leayis to a number of geoaetricaiiy 
realizable physical invariants of the orbxts^ Consequently, aost ui 
the theorems rely only on simple geometrical relatxon^hips* Despite 
Its siiplici.ty, this planetary geometry is powerxui enoui^jn to treat < 
wide range of perturbations with relative ease* It xs leit' that this 
"rreataent provides a better view of "what doing physxcs reaily 
like" than the standard roate via algebraic aanipuilatx<3hs* Ihe 
doca»ent concludes with suggestions ^or further research xnto tne 
qeoffietry of planetary orbits. (MP) 
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develop a theory of orbits ^ ^ .e^se-sr.jare zentral force law which 

differs considerably fron^, the ^zja -..j >e a3::ro5,,h In particular, we make 
no ej<plKit'%se of calculus by ber^^^^5 ^^t^ qi.v:^itoMve aspects of solutions, 
we are led to number of geor-etn ' ; ^p^^zcibV pnys cal invariants of the 
Orbits. Consequently most of o';r t^i-.r^'-s ^--V on'y on sinipie georretrical 
relationships, ^- - — ^ 
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enough to treat a w^de^^ange of ^ert.'' 
without introducing any rore ma^hu'^^y 
ihe^paper concludes with suggestion: v 
of planetary orbits 



.^ne'iry ^^^^onetry is powerful 
'lO'^'"^ reldOwe ease Furtherr-ore, 

'nta-r ^^i; ]uantitative results 
^^j^^ne: "^esear h I'^to the geometry 
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h Some Qualitahve KesuUs 



1 Intfoductton 



if/c 

f ■ » 

From junior high sch^ on, students of science are taught that Kepler's Laws 
describe the motion of planets around the sun They are given no hmt of how they 

V 

'themselves can 4jnderstand the "why" of these laws By high school the stjdenls have be?n 
taught Newton's discovery, that the inverse squarei force law accounts 'or those beautiful 
ellipses, but the connection is not yet for their eyes After a year or so of college it s finally 
time to plow through the thoroughly siandardued and unmotivalea proofs, using intricate 
manipulations with differential equations • 

In this paper we outline an approach to orbital mechinics which is accessible to 
beginning physics students and presupposes no knowledge of calculus We give an 
elementary (yet mathematically correct) treatment of Kepler's Laws and also investigate a 
simple first order perturbation theory for orbits in an mmse square field Our theorems 
and proofs arise natiirally from irymg to understand orbits in terms of their physical 
invariants therefore ^eel that our treatment provides a better view of *^hat doing 
physics isVeally like" than does the standard route v.a algebraic manipulations 

^ The ktV to thf method lies in consdenng the velocity space picture for a planet's 
motion about the sun The concept of a »«locity space is not normally encountered by the 
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uju^lly c^nu:ie''*J '^a iJv^ri'.c .*[ '■\,*, sp'^^t pi: urc i^'^rt/cle's motion 

h^s rh*^ co'^ ^'r.niin rn* h^'-' 'h^: qualrvav* content of Nfwton$ F-m?^ 

tor our p'i^|jc »- nA^-FA'' i> sii^'pU T'*"' i'^' m^ff^'tinn. hnwten objects t^k^platce 
by ^ no:Mh^^:ii^n c y^^'lcrrt ^^'H*'- ^'^^n hr -i ch*^n?f in po I'Eon Appreciation of thi$ 

ip^Cf^,is H rn u"^! 'rn] * or ^ , ^jirjr-in^ r^' on ^ pNfnaril r f oncf pp;^ I breaMhrough In Pari 
HI when wf '^ji'.uT^^'u'b^Mon we sh^ll wh-^' nrh di^idtncis cun be y^ielded by 
looking at fjhvsr^l phpncwn^^n^ in ^1-^ ^-r^h' cGnc^p'u'^l fran? in this case, velocity space 



In ^h<^ rollu^irir fjr*^>* n'^ \< at * jfd to walk rather a narrow path between 
two pxtr^rnvs Cjf\ one h^ud a d^^urjfjMon of our rntrhods and results would take no more 
^han a p"^?,^^ \' us^^d ihf full prtosjQn of nathematKal apparatus (mcludmg, 
ralcijlu^i ^vaJl3b)p ro <vM*nff sturlrnf^ r ^ t<^/^ yf:ars of univer;[^ty education On the 
othpr h;^n'l wc > t)ij|rl hv.f ^p»ni rormjtrdbly morf space d^^velopmg a complete and self - 
ron^^jnfd course for^very fatjy phy^ici s*udKnts> Since wf^ feel the niaterial can be^useful at 



boih Ip^f'K uf phyULi^ tuULaUun. we have aiienripied a compromisf We apbiogize both to 
ihose who Imd oui preipn^auon eMrnd*»d and peihapj verbose, and to thoic who might 
find It sketchy and incomplete 
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Wc gratefully acknowledge the inspirauon and Bncouragement of Seymour 
Papert He introduced us to this way of thinking aboul orbits, and pointed out the basic 
results described in Sections 2 and 5 These sections closely follow parts of his paper 
(reference 1) which presents a broader view of the conception of education in science and 
mathematics from which this work grew We would like lo thank Ihe editor and referee 
from Tht AmerUan Journal of Physics for many encouraging and helpful comments and 
also to thank Suzm Jabarl of the M IT Artificial Intelligence Laboratory for preparing the 
illustrations for this paper 

2 The Orbtl is Closed 

/ 

, Standard approaches begin wit)i the arduous task of proving that planetary 
, orbiLs arc precise ellipses *We begin by proving a more qualitative proposition, that orbiU 
arc closed In doing so we dispense with a great deal of analytic clutter, and the important 
special nature of inverse-square orbits which makes them closed comes into central focus 

We will prove that no orbit like that In Figure 1 is possible 
> Jf a planet crosses a half-line from ihe sun iwlce, ihen II crosses ll al 
Ihe same poinl each ijie-nol furiher oul or closer In 



We assume two pieces of knowledge 
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1 The force on the planet, when it Js distance r from the surf, Is K/r^ towards the 



-sun 



/ 

2 Angular momentum Is conserved We use this in the form of Kepler's Law 
that the radius from the sun Co a planet sweeps out equal areas In equal times This can be 
easily derived and we remind readers of iis simple geometric proof in the appendix 

f 

Now consider diametrically opposite pieces of the orbit which subtend the same 
(small) angle measured from the sun, as in Figure 2 Kepler's law states that 



A2 ^ta 

What else do we know about the area or tt^e time? Geometry tells us that 
. 2 



r 



Those r^'s are too suggestive for us not to make a connection with 
F - K/r^ In fact 



2 

hence — ! — - L 



S . . ^2 A, zit 



an<J we conclude F|Al| - Yf^LK^ Since Fj and F2 pull in opposite directions, 
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We can id^nttfy these terms Accoiding to Ne^'on's Gecond Law,^ on each piece of 'ihp 
orbit, FAt IS precisely Av the change in velocity^ whtch we call the "^kic^ associated to 
that piece Thus the last f^quation thai: che change fn velocity over one piece of orbit 
exactly cancels the chinge on the opposite piece Starting on the half-line which the orbit 
crosses twjce, divide the orbit a!! the way around into similar pairs of opposite pieces The 
total change in the planetN velocity bf^ween successive crossings of the h^lf-line n the sum 
of th^ ch-m^ts in velocity ov^r each small piece, adding t>iese up in opposing pairs, we see 
t^^si^the to^>l ch inge is z^ro Wh^n^ver the p'.n^r rrossej a given halfHme, it has the same 
velocity ^ 

m 

Now Kepler*s dictum of equal area tn equal time allows us to conclude that at two 

crossing?: of the half-line, not only velocity ^3ut distance from the sun is the same Figure 3 

Shows the areas 6wepC out by the planet ii\ some short time At after Successive crossings of 

the half line through A. B, and O The velocities at A and B are equal Therefore the 

pieces of orbit AC and BD are both equal totAt, but the area of AOC must equal the 

area of BOD Then A equals B and the orbit closes 
« 

3 A Theoren tn Velocity Space 

The preceding proof rested on the fact thit the kicks over opposite plecei of the 
orbit have equal magnitudes 



Ffgure 1' AC=8D, area AOC=area BOD, therefore OA=OB 57T3^A=B. 
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But, to derive this equation it isn't neces^^ary th?^t tht pieces of orbit be opposite We need 

only ^'iX j 

^ 2 

A, 

— = — T '> 

^ ?nd thii. is'^irue for any two pie<.e\ of the oibit over which the radial Jungle changes by the 
same smalt amount So, if we diviae liie v.^o^*^ oibit into small pieces subtending the same 
angle , the ''kick" vectors for the various pieces all have the same length (Figure 4) Not 
only are aij the lengths equal, but the lotation between sucessue kick vectors is constant ajnd 
equals Atf 

Tfius we have a very simple algorithm for generating the changing velocity as 
the planet moves along its 01 bit Starting at a given velocity vector we add on kick vectors 
one ^fter another Each addition is a step of constant length and successive steps differ in 
direction Ijy the constant turn. Ali It Is easy to see that the algorithm GO FORWARD a 
short distance, TURN through a, small angle, CO FORWARD the same short distance, 
T9RN through the same small angle, , will generate a circle We conclude that the Kick 
vectors line up along a circle (Figure 5) 

Wc can interpiPt "adding on succrssive kick vectors" by introducing the notion of 
vd^fcity ipate The velocity^of an ubjpct i> usually described by a ""velocity vector*, th^t is. a 



r 

dtrection and a sp-^ed (Itngth) In comparm^ diffcrtrnt velocities it is useful to put the Uil 
of all velocity vectors down at ^on^ comman^.point, 0, and to depict a velocity by the point 
where the tip of the velocity vecter hnd$ jwuh this convention, we can draw two different 

4 

0 

pictures to describe the motion of an object 

(I; the colit^cuon of successive positions of the object m "reaf space, and 
(2) the collection of successive positions of the tip of the object's velcxity vector 
'path in "velocity space", ^picture of how velocity changes 

The second picture js called the "velocity space path" or v'velocity diagram*. 
Figure 6 exemplifies these two kinds of diagrams 

Velocity IS the thing that chanrges position, kicks arc the things that change 
velocity To get from an object's position at one instant, t, to \is position at t*At, we add on 
the Sector vA[ To get/rom t^e object's velocity at one instant to the velocilylit a slightly 
later time, we add on jB^^^^^*^^- Adding up successive vAt vectors gives the 
position space path; adding kick vectors gj^^es the velocity space path We can now. restate 
our result as 

CircU Theorem 

For an cbject moving tn an tnmse-sguare field, the vitocitj ipaci 
path liei on a circle 
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To avoid confuston w^^ point out that the center of this Circk i$ not necessarily at 



the.origm in velocity spa^ 



4 The Velocity Space Path 



Our Circle Theorem' 1*;^] Is us that the v^locif^ space path heson a circle But i$ it 
a complete circle or just pari of \(> We can answer thai question and a bit more 



showed thai an orbit which does manage lo get all the way around the sun 
presses every half lme from the sun exSictly once Such an orbit is a simple closed curve 
In a complete i evolution. iberero:e. the direction of the planet's velocity vector must change 
through 360* (See refer^ce 4) That means thai in velocity space ^oMht path meeu every 
half-line from the velocitj^ space origin Ii follows that closed Orbits in position space 
correspond to complete circles in vefocity space, and we have learned, besides, that the 
origin of velocity space f^s irjygde the circle 

We now have a good qualitative picture of the velocity space path for a closed 
orbit (Open orbiu are diKus^sed m section II ) In Pait 11* we will extract Information about 
the position space orbit from our velocity space diagram 



I 



f 



PAGE 12 



II. Invananfs ol ihz Qvhii 



5 Angular Momentum 

4 

As we have s^en/the velocity space path of planet in an Inverse-square field hes 
on a circle One obvious invariant of a circle is its radius How can wc Interpret this 
invariant physically^ 

We got the circle in section 3 as the result of the algorithm ^forward distance D> 
turn angle Aff. repeat" As one can see {Figure 7), this generates a circle of radius In 
our case was an arbitrary small angle and D was the magnitude of the Kick F^t over 
the corresponding small piece of>'orbiU LeLung u denote the radius of the velocity circle, we 
have ' _^ 



It u not immediately obvious that this is a constant However, we can simplify the 
expression using the fau from geometry ihaL the area iwept out over a small piece oT orbit 

IS 

A-ir^Ae 

2 

Then u-E4i" Fr^At We can eliminate the apparent dependence on the non-constant term 
by ustpg F-K/r^, we obtain % 



u 
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V 

Th«'term 2A/At. which telU how fast the planet is sweeping, jout area, is precisely the 
cbnstant called angular nomentum, L (See Appendix ) Therefore 

The radius of the velocity circle equals the force constant 
K dividtd by the angular momentum L u-K/L 
For a fixed gravitational field, the radius of the velocity circle tells us the planet's angular 

momentum a larger radius gives a smaller angular momentum. 

/ 

6 Ortentatton 

The velocity circle has another invariant so obvious it Is easy to overiook-the 

position of lis cenler in velocity space As we remarked above, ever, .hough the origin 0 Is 

tnslde the velocity circle It need not be<i/ the ctr.ur of thr tiVde Let I be the vector 
♦ 

running from the origin in velocity space tp the center of the velocity circle, and let u be 
a radial vector of the circle (Figure 8) In terms of i and u. we can think of the planet's 
path in velo*.ity space as follows at each moment the velocity v Is the sum of a ccnuant 
Victor I and a vector u of constant. Ungih (equal to u-K/L), v-z*u The velocity space 
path IS generated as the radius u sweeps around the tip of the Invariant vector^ 

There is a quite remarkable relation between the motion of f. the position space 
radius vector (tall at the 5un, head at the planet), and the motion of this "velocity space . 
radius.' u 

Correlation of Angles' tn Position and Velocity Space ' 
« 

At each homent the planet's radial vector r f$ peipendlcidar to the 
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radtus u of th^ velocity circle 
To see this wc examine how Ihe kicks fil into both diagrams In position space each kick is 
parallel to the radial vector r In veloaly space the same kick is tangent to the velocity 
circleand hence perpendicular to the velocity circle radtus u- 

Hence-u is perpendicular to r 

Correlation of Angles is a powerful principle It tells us 

(1) Each point on the velocity space path corresponds to a unique point in the 
planet's orbu (The planet cannot attain the same velocity at two different points in its 
orbit) 

IT 

(2) The planet's position vector sweeps around the sun at the same rate and with 

th^ same direction (clockwise or counter- clockwise) as its velocity sweeps around the circle 

■ 

in velocity space The two are always 90 degrees out of phase (Figure 9) 

. / Now we 'can, give more meaning to the i vector The planet's speed, v, is 

r 

* greatest when the u and i vectors are lined up. least when they arc opposed v^^^-u-z, 
v^^^^mUii Therefore the z vector points in the direction of maJ^imum speed and 
opposite *to the direction of minimum speed 

The points of greatest and least speed occur where the velocity vector is parallel 

■ 

to u.arhd so perpendicular to the position spa^ radius r It Is not hard to show then 
thai the point wheie speed auains ils ma?;im*>m (respectively, its minimum) corresponds to 
the minimum (respectively, maximum) distance from the sun The reader can fill in the 
^ details of the proof sketched below 

*- A 



POSITION 



VELOCITY 




POSITION 



VELOCITY 
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POSITION 




VELOCITY 



Figure 9 Snapshots of position and velocity over an orbit Ujf. 
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Step 1 At % maximum or minimum distance the velocity 'w must be 

perpendicular to the radius r 

Step 2 The velocity diagram shows thatthere are precisely two points where this 

can occur \ 

^ Step 3 Conservation of Angula^ Momentum Implies that r^^^j^ corresponds to 

^min r^^, to v^^^ 

Now we have some more qualitative information about the shape of the orbtL 

there is preciseljf^one point of maximum distance from the sun, and one of minimum 

distance They occur on opposite sides of Ihe sun since the corresponding^ vectors point 

in opposite directions The r vector determines the orientation of the orbit. It points In 

<* 

the direction of maximum speed (FHgure 10) 



7 Zhapi and Symmetry 

What more does the length of "r tell us about the orbit? Consider what would 
happen if f vanished There would then be no direction picked out for maximum speed 
or distance from the sun The planet would have to travel around the sun In a circle at ' 
uniform speed (Another way to see this v would be equal to iT and therefore always 
perpendicular lo'r and of constant length the characteristic of uniform circular motion ) 

This suggests that i indicates how the orbit deviates from a circle. We can 
make this precise One obvious measure of the non-circularlty of the orbit Is the difference 
in the extreme distances from the sun, 

' > " 
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If we want an invariant that depends only on the shape and not the size of the orbit it is 
better to see how much (he ratio ' j 



'^max^'^mln 



differs from I The length of i, relates the maximum and minimum speeds 

^mln-"-^ ^max-"*^ 

To relate speed to distance from the sun we use angular momentum. If the angular 
momentum is L then 

^*^min ^ '^max^'^max ^ ^min 
because at these places in tjjje orbit v is perpendicular to r (Section 6) Therefore. 

Tmax ^ Vmax . ^-^^ 

Since u-K/L we can also write 

*"max Kt Lz »> 



But K depends only on the nature of the gravitational field so we see that our "shape 
invariant* is determined by Lz The larger Li, the more the orbi( deviates from a circle 
In section 11 we will derive the analytic result thai the orbit Is an elUpse and Lz determines 
Its eccentrtctty 

In fa^r knowledge of K and the 'velocity diagram essentially determines the 
planet's motion The maximum velocily can be read off the diagram immediately, as can 
, u-K/L, so we know L The shortest radius in position space has length ^rnin'^'^max* 
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IS perpendicular to v^^^^ (Section 6) Having this one vector, r^^^, vre can generate the 
orbit starting at the positiOtn determined by t^^^ with the following algorithm 

1 Travel a short distance vAt in the direction of v 

2 Measure the change in angle, Aff, in position space, 

3 Find the velocity arthis new angle (by rotating u through Afl and consulting 
the velocity diagram) 

4 Return to step I This generates the entire position space path 

Notice that the velocity diagram is .symmetric about the Ime determined by ^. 
The above algorithm translates this fact into a symmetry of the position space orbit 
Starling at the nearest point to the sun, construct the orbit in the forward direction for a 
while, along 'vj for Aflj. then along "vj^ for Afl2. and so on Now go back to the 
starting poini and run the algorithm backwards with the same sequence of Afl s Since the 
velocity diagram Is symmetric wc generate \\\t iamc small segmcrUs of orbll, except that they 
have been, flipped about the line perpendicular ioi Therefore the entire orbit is 



symmetric abouit this line ^ 

\ 



8 Summary 

We have so far obtained the following information from the velocity diagram. 

1 The radius of the velocity circle detel^mlncs the orbll^s angular momentum- 
u-K/L \ 

2 The center of the velocHy circle detcrmlnej the orientation of the orbit (T 
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potnti in the direction of maximum ipeed) and It* "shape" (Li determinei rmax'''mln^- 

3. From the velocity diagram, we can aIgorlthm)ca!ly determine the whole orbit. 

V 
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9 Tht Ptrturbatton Formula. Radial Thrust 

Ic IS In the Jtudy of perturbationj. or how orbiu change under small kicks other 

lhan those given by the sun, thai our use of velocity diagrams really pays off. It pays off 

for a very good reason, which we mentioned in the introduction as a qualitative form of 

t 

Newton's Second Law of Motion ^ 
Force acts on the paths of particles by changing velocity and not 
posiiion 

If we fail to take account of this fact we may be faced with situations tha( appear 

counterintuitive For example, suppose a spaceship In a circular orbit around a planet 

S 

apphes a small outward thrust (Figure ]la)^ 

How will the orbit change? "Intuition' may suggest that the orbit will elongate in the 
direction of the thrust^ something like Figure lib In fact, the orbit will elongate, but in a 
direction ptrpendtcular to the kick as in Figure lie 

^ To Ainderstand this we consider how the kick changes the velocity diagram The 

spaceship started in a circular orbit whose velocity diagram is centered on the origin 
Since force affects velocity and not position, ll is reasonable (and we shall show below) that 
the effect of the kick in velocity space really ti just to move the velocity circle in the 
direction of the kick (Figure 12) The corresponding change in the position space orbit is 
the "counler-lnluitive" effect described above 

Our strategy for studying perturbations will be to see how kicks change the 





^ " Figure 12 The perturbation induced by an outward kick, 
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velocity diagram More precisely, we know lhat the shape and orientation of the orbit is 
determined by Lz, the z vector times the angular momentum, so we want to fmd the change 
in Li, A(Lz), produced by an arbitrary kick 

* The banic velocity space equation, v - ? ♦ u, gives Lt L v - Lu Since u ha$ 
length K/L we have Lu ' Ks) wheret is a vector of unit length whose direction Is 
determined solely by the object's position (It is perpendkula/to the radial vector^) To 
compute the effect of a kick on Li -i L v - Kt, notice that since kicks do not affect 
position, t IS unchanged K is aUo unchanged Therefore the changein Lz is the same a$ 
the change m Lv» and the first-order approximation to the change in a product of 
changing quantities gives 

Ftrturhatton Formula A(Li) -"^AL ♦ LAv^ 

/ I 

We can use the Perturbation Formuta to tidy up our discussion of the *'radlal 
thrust problem"^(Figure H) Since Av is in the radial direction, the angular momentum 
does not change (AL - 0), so the formula implies A(Lz) - LAv. This means that the 
velocity diagram changes from at vector of zcw to a*! vector in the ^ircctlon of Av 
(Figure 12) 

Intuitively, the velocity circle is "pushed" in the direction of the kick Note that an inward 
kick at the bottom of the position space orbit wouid have the same effect 
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\0 Thngenttal Thrust, Solar Wind, i^^*^^ Lam 

X 

«% 

\ 

In this section we^ppty the Perturbation Formula to some other orbit problems 
Tangenttal Thruu Suppose again that a rocket starts in a circular orbit, but this 
time firovides a tangential kick (Figure 13) To determine A(Li) - "v^L ♦ LAv we note that 
vAL is an impulse in the direction of A'v since Al is positive and "v Is paratW to Av 
Then the newly created z vector must be in the direction of the impulse (Figure 13) The 
elongation in position space is again perpendicular to the kick 

The Soiar Wtnd'^ Assume th^ rocket is affected not only by the planet's gravity 
but also by a small constant force (Luehrmann's "solar wind*) If the perturbing force is 
smalt compared to gravity^ each revolution of the rocket will be nearly an ellipse. We can 
therefore think of the orbit as an ellipse Jghich varies througj^time To compute how the 
ellipse changes we vie,v the wind as providing impulses all along the orbit (Figure 14) and 
sum over one revolution 

( The UAt contribution is a net change In the direction of A^ To compute "vAL 
we notice that AL is fjosuive on Ihe bottom half of the orbit and negative on the top half 
as $ho\^ m Figure 15b We can sum the"vAL's by exploiting the symmetry of the orbit 
The vertical components of the'vAL's on the left cancel the vertical components of the 
"vALs on the right, leaving only a horizontal component In the direction of A'v' (Figure 
I5c) This adds with LAv to produce at' vector In the direction of the wind Intuitively, 
the velocity circle gets "blown" In the direction of the wind The orbit elongates 
perpendicular to the wind as in Figure 16 




t 



V 




c 




The vectors v (a) and v/iL (b) fqr the solar wind. . In (c) 
we see that the vertical components of l^'AL on the left 
cancel the vertical components on the right 
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Since the wind does not affect the symmetry of the orbit about the vertical axis we can 
apply the same analysis as abov* to show that the iTz vector continues increasing in the 
direction of the wind From this we conclude that the orbit becomes more and n^ore 
eccentric while the direction of t remams constant^ and the orbit continues to elongate 
ptrptndUular to the wind 

The orbit becomes closer and closer to a straight line, and it eventually reaches a 
point where the "smalf wind can hav^rge qualitative effects over a timescale of ies5 than 
one revolution (the orbil in fact reverses directionX and our method of averaging over an 
entire revolution becomes inappropriate 

The r'^^**^ Forct Pjitld If € is a small constant (we will take it to be positive), we 
can treat the central force field of magnitude r"^^**^ as a perturbation of the r'^ field The 
perturbing force is some force (positive or tiegative) In the radial direction To understand 
how this perturbation affects the orbit, we maJ^e the important observation that the shape 
of a r"^ or r"^^*^^ orbit (3oes not depend on the scale which we use to measure radius 
Therefore we can determine shape by using any scale which makes it rnnvenlent to 
compute the effect of the pf iturb-^lKiji I-oi the orbit shown in Figure I7» we scale to make 
the distance OP equal to one Since i/r^ < 1/r^*^*^^ for r < I and l/r^ > 1/r^^*^) for r > I the 
perturbing foicc is as shown 

For this perturbation the kicks are radial, so L is constant This means A(Li) - 
LAz! but from the perturbation formula, AL ' 0 implies A(Lz) LA^ Hence A^^- At* 
Now wf can sum ta over an entire orbit The left right symmetry of the orbit and 
perturbing force m^ans that the sum of horizontal components of the kicks must canceK 

3 ; 
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and the f^l At is downward, thai peipendicuhr to the original t SubsequeJ^Ats will be 
perpend wbf to the current f and this results piimanly in a rotation of ^, not a change In 
^ iength (Figure 18) The "rf>ajor axis' of the position space orbit, by the consequeaces we 
derived fror^orrelauon of Anglermust follow this counter-clockwiseN^otation Though 
the orbit retails its shape hi. U preces.^es (Figure 19) ^ 

1^ Warning It should be rem9irked that in%e preceding two examples we looked at 
the F vectors as representing Av for the perturbation formula Of course we should have 
used TAt. but. because wt the symmetry involved, the At factor can be ignored in those t-ujp 
cases In more complicated situations, though, this does beconr>e an issue For example, we 
invite the reader to use the techniques yjt this itUion to tieat the p rturbalion mduced by 
an oblate^fiun 



: II The Orba is a Conic Section 



An objection that is sure to occur to some of our readers goes something like this 
"All these mtuitive methods are fjne. bi/t if you want useful quantitative information you 
have to return to the standard differential equations you've been trying to get along 
withOLit" ^ 
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Of course there are orbit problems our simple methods won't handle As for the 
standard results, however, we are able to derive the orbilal equation directJy from our 

r 



velocity diagram^sing no more than trigonometry. 



The orbit is describe in polar coordinates by the equation 

u *-2 cos^ 

The proof is a natural correlation of the basic quantities, rl'v, u , 1 and L using 
the definition of angular momentum At a point In the orbil when theU vector and thc"z 
vector differ in direction by an angle $ we construct the angular momer^tum triangle (sec 
appendix) 



The area of the triangle OAB in Figure 20 is t^defmitlon L/2 If h is the height 
of the triangle tH^n * 

2 2 

Since u and r are perpendicular, the height of the triangle is given by h - u - z co$ 5. 
Therefore 

I* L - rh - r (u ' I cos ^) 

Here B represents the angle in ueloiuy space between ^ and ihc fixed vector l Correlation 
of Angles implies that 6 also measures the angle In position xpace from to a fiXed /cctor 



POSITION 




VELOCITY 



POSITION ^ 

Figure 21' Sample ^bits in position and velocity space. 
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perpendicular to 2 Therefore r and 6 are polar coordinates in position space 

The above equation describes a conic section When the origin of velocity space 
lies within the circle, u>z and the orbit 15 an ellipse When the origin is outside the circle, u 
< 2 and the orbit is hyperbolic When the circle passes through the origin, u - 2 and the 
orbit IS parabolic (Figure 21) 



Writing the equation ^n the form 



2 



V 



r 





w« get the "standard form" for a conic section and see that / ^ the eccentricity, and 
L^/K IS the radius of tKe orbit when the eccectricuy is zero 



12 Co7is(Tva:toyi of Err^rgy 



Energy conservdfion dops not a^ise naturally using this geometric apprcB^ch 



although we can obtain the result as a simple application 



Apply the law of cosines to th^ velocay diagram in Figure 22 to gel 



1.. 



Figure 22: = + - 2u2 cos^?. 
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2 _ 2 2 " 

V "2 + u -2uz COS d 



SStitutlhs: z coxff - u - y (f rom Sectior^ 11) we obtain 



2 _ 2 2 2uL 

V - Z -LT^-^^ 

r 



and henc*; (since u * K/L) 



1 

V ^ z'— u 



K 



2 r 



Since I and u are constant for the orbit, 2 ^ constant, llr^ total energy. E It is 

interesting to note that when the planet crosses the seml-nlno*^ axii fTperptndicutar lo"^ Tf. 

vV 

T and V form a right triangle wuh v^-u^ i.^ hence the kinetic energy /2 1* exactly the 
negative of the tocal energy there 



J3 Kepltr's Third Law 

We can use the relation of angular momentum to area swept out. 2A-Lt. to 
compute the period of the planetV revolution In one period the planet sweeps out the 
entire area of ijj/cfliptical orbit The are^i of an ellipse of semi major axis a and 
eccentricity e is given by A - aV^-e^ For the orbit we have 

1. 
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a - 2 ( "^min + 'max) = - ( + ^ \ 

s\u — 2 u — z } 

= ^" - K 

a 2 
U - 2 

The eccentricity IS e - i/u, JO \J 1 -e^ = ^ = 



^ 2 ^ ^ 



Then the period is d«?termmecl by 



u 



Of 



L u 



In terns of qu^nti'ies appearing in the ve'oci'y diagram we get 

V - 



■ 

Rath#>r than tre^t the hyperbolic case in detail, we leave the reader to verify the 

following 

I For an open orbit, the arc of the velocity arcle which is actually traversed is the 
part shown below, bounded by the tangenis lo the circle through the origirt of velocity 
sp^ce Velocity space geometry gives the correlation of energy with limiting velocity, v 



1 . 



•J?^ JlT (Figure 23) 



2 The deflection angle (angle between the two ajyrripiotej of the hyperbola) can 
be. easily found as a funcUon of energy and angular momentum tan m\mJ^K^ (Figure 
24) • 



15 Su^^emorii for Further Research 



We have by no means exhausted the study of the jgeometry of orbits In this 
paper The geometry of orbKs, particularly the perturbation theory, 1$ a rich source of 
problems, even of mini-research projects of ihe lype described by Luehrroann . Below we 
maVe some suggestions for problems and study topics 

An tnstruUioe paradox The Galilean uansformalion requires the rcjalion between velocities 
measured in twQ different frames of refrence moving with relative velocity Vq to be v - v 

♦ Vq This It IS a Simple matter to mo/e into d frame with relative velocity z and transform 

/ 

the"z vector for an elliptical orbil to 0 Why then does one not observe a circular orbit in 
the new frame? In particular, whal fails in the algoiilhm uf SecUon 8 which doet generate 
a circle In position space given a circle centered about the origin in velocity space? 
An Aid to Astrogatton? Suppose we had to pilot a spaceship in a gravitational field (such as 
simulated in computer 'space war" games) Would a velcKity diagram be a useful addition 
to our instrument panei^ Foi example, to change from an elliptical orbit to a circular orbit 
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■ 

■ 

r 

■ 

we need only cornult the vetociiy diagram and apply a force to cancel the z vector. On the 
other hand, lots of information is lacking if we use oniy the velocity picture Intercepting 
another spaceship is a tricky problenfi involving timing (Although merely matching its 
orbit IS easy) What other Instruments should supplement or possibly replace a velocity 
diagram^ 

Geometry of the harmonu osullator A lundamenial geometric property of solutions to the 
Ji/^ differential equation is that they have a vector constant of motion (zior the maximum 
velocity, or the Runge-Lenz vector are all[ possible choices for this constant) The two 
dimensional harmonic oscillator wHh equal mode frequencies has a similar structure. 
Solutions have an obvioui axis which may alio be awigned a magnitude in any number of 
ways Can one develop a useful velocity space geometry for that system? Can one treat 
simple peVurbalions, as with orbits? 

More Soiar Wind A further discussion of the solar wind phenomenon could make use of 
the fact that the force field is conservative and therefore the energy of the orbit Is constant 
This implies that the angular momentum decreases as the orbit becomes more eccentric 
Thus the velocity circle is not only 'blown by the wind** but also the radius becomcs^larger 
and finally infinite when the orbit degenerates to a line Show that th<; changing velocity 
circle always passes through two fixed poinu in the plane (Figure 25) 

Our method of averaging over entire orbits Is only a first-order perturbation 
theory whereas the fornlula 
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(s IS any parameicr) is exact One would hKe to have a treatment of the solar wind which 
^ works near the turnaround points, zlu approaching I 

Finally, is there a complete perturbation theory based on the geometry of orblli? 
In particular how can one treat perturbations out of the plane of the orbit? 
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Appendix: ;^n^u/dr Momentum and Ktpleri Second Law 

/ * - 

ThroMghout this paper we fiave been assuming Kepler*s Second Law There Is a 
A Simple get)mecrl^ prwf of tjbis which we place here In an appendix because u did not 

; ^ ^ 

originate with us it can be found In Newton's Principla^ 

^ . ^ ' ^ 

The atigular momentum, which we denote by U is defined to bt twice the area of 

4 

iang^le determined by^he velocuy vector and the radius vector from the sun to the 

planet As shown below L is constant if the velocity doesn*i change - the triangles have 

^ual ureas since they Ijave equal bases (the length of "v) and equal height! (FigUte Al). 

Morf remarkably.'X remains constant tf we ci^ange the veloctly bj applying any kick in the 

radial directton (towards or away from (he.S}in) The effect of a kTck Av on the angular 

momentum triangle is illustrated bcloAi (Figure A2) The kick changes'? to"v^ but triangte 

Operand OPOj have th? same base, OP, and the same height (h in the dlagramj^smce OP 

and QjO; aie parallel Therefore OPQ^and OPQ; have the same area, and angular 

momentum is unchanged * 

^ A planet moving about the sun, subjfct to no force but the sun's gravitation, has 

every Vick applied in the radial direction None of ihese change the angular momentum, 

/ 

. whrch^is therefore an invaiianl of the pjanel's oibit To find a geometric interpretation of 

' H • ^ ^ A 

this fact, m e«tamine the orbit at time intervals Al small enough that the velocity does not 

change much over each interval In each interval ibe radius vector sweeps out a ^mall 
« 

triangle . The area of one of these smjll triangles u ' >- 



• 



0 



o 



Figure A^; Area OPQ = Area OPQ* 



J 



1 ^ 

< I 



'^h ai iriuPig"''^ ^^^^ ^^^P' ^^""^ '""^ 



T -At, ♦ ♦ 

IS ths ium of the areas of the small trianglei 



I 1 



= l-LT 



This gives Kepler's Sc<ond Law ^ 
For a body moving m a radial force field, the radUl verfor sw^^s 
oui equal areai in equal time>. 
It n unfortunate that this proof is not more often presented Jn phji^sics coursei (although 
Feynman ^ discusses u and there is a movie ' demonstrating this argument) 
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